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ABSTRACT 

We measure the physical properties of a local multi-component absorption-hne system at V0 
~200 km s"i toward the quasar PKS 0312-770 behind the Magellanic Bridge (MB) using Hubble 
Space Telescope STIS spectroscopy in conjunction with photoionization modeling. At an impact param- 
eter of ~10 kpc from the Small Magellanic Cloud (SMC), this sightline provides a unique opportunity 
to probe the chemical properties and ionization structure in a nearby absorption line system with a col- 
umn density of log A'h i ~20.2, at the transition between Damped Lyman Alpha (DLA) and sub-DLA 
systems. We find that metallicity of —1.0 < \og{Z / Zq) < —0.5 and ionization parameter of —6 < \ogU 
< —5 for three low-ionization components and logU ~ —2.6 for one high-ionization component. One 
component at Vq = 207 km s^^ shows an a-element abundance log(Si/H) '^—5.0, making it ~0.2 dex 
more metal rich than both SMC H 11 regions and stars within the MB and the SMC. The N/Si ratio 
in this component is log(N/Si) = —0.3 ± 0.1, making it comparable to other N-poor dwarf galaxies and 
■~0.2 dex lower than H 11 regions in the SMC. Another component at Vq = 236 km shows a similar 
Si/H ratio but has log(N/Si) — —1.0 ± 0.2, indicating a nitrogen deficiency comparable to that seen in 
the most N-poor DLA systems. These differences imply different chemical enrichment histories between 
components along the same sightline. Our results suggest that, if these absorbers are representative 
some fraction of DLA systems, then 1) DLA systems along single sight-lines do not necessarily represent 
the global properties of the absorbing cloud, and b) the chemical composition within a given DLA cloud 
may be inhomogeneous. 

Subject headings: Magellanic Clouds - galaxies: individual (Large Magellanic Cloud, Small Magellanic 
Cloud) - galaxies: abundances — absorption lines 



1. INTRODUCTION 

A physical connection between the Large Magellanic 
Cloud (LMC) and the SmaU Magellani c Cloud (SMC ) 
was implied as early as 60 years ago (jShaplevI |1940D . 
A continuous H i gas structure between the LMC and 
the SMC, know n as the Magellanic B ridge (MB), was 
first reported by iHindman et aP (|l963f ). ar id it is found 
to be at a distance of ~50-60 kpc (e.g., iHarries et al.l 
120031) . There is substantial evidence for a young stel- 
lar population in the MB, with s tars as young as ~ 
20 Mvrs (Jrwin et all 119851 Il990l: iGrondin et all Il992t 



iHamblv et all 1X9941: iDemers fc Battinellil 1199811 and with 
mole c ular clouds detect ed (jLehned 120021 : iMuller et al.l 
l2003t iMizuno et al.ll2006h . These young populations were 
probably formed locally, because they are not old enough 
to have escaped from the next nearest star-forming re- 
gion (i.e., the SMC) based on their peculiar motions. 
On the other hand, the larger-scale tidal structure that 
is connected to the MB, the Magellanic Stream, con- 
tains only gas, with no evidence of a stellar popula- 
tion (jCuhathakurta fc Reitzell [1998^) . The physical con- 
nection of these structures are still un der discussion 
(jKallivavalil et al.ll2006l: iPiatek et ai]|2008D . 



Although the formation mechanism of the MB is not 
understood in detail, gravitational tidal inter actions prob- 
ably played a primary role (e.g., Gardiner et all 11994 
iMurai fc Fuiimotolll980f ). Numerical simulations indicate 
that the MB was created through a close interaction be- 
tween the LMC and the SMC about 0.2 Gyr ago (e.g., 
iGardiner fc Noguchil [1991 . Thus, the MB is an ideal 
target for studies of the influence of dynamical interac- 
tions on star-forming activity. The MB also provides a 
unique opportunity to study nearby star formation in a 
low metallicity environment. The MB stars have metallic- 
ities of only one-tenth that of normal Population I Galactic 
stars (i.e., even smal ler than the SMC stars by ^^0.5 dex; 
iRolleston et al.lll999f) . 

The abundance pattern of the MB has been stud- 
ied through stellar populations and gas clouds that are 
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Lehneil 120021: iMallouris et all 120011: IWeltv et al boOl 
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cently, Carrera ct al. (200^ found a possible metallicity 
gradient in the SMC toward outer regions from ~1 to ~4 
degree from its center. However, stars and circumstellar 
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regions could be biased because their physical conditions 
are significantly influenced by nearby stars. Analysis of in- 
terstellar gas should yield more representative properties 
of the MB and its abundance pattern. 

Quasar absorption lines are powerful tools to investi- 
gate the physical conditions of absorbers located along 
sight-lines to quasars. These absorbers are detected in 
quasar spectra, regardless of the source luminosity, which 
enables us to collect a homogeneous sample of absorb- 
ing clouds from early cosmic epochs to the present. Al- 
though there have been several attempts to trace the MB 
gas along the sight-lines toward quasars, these were based 
on radio-spectra that cover ed H i 21-cm absorption (e.g., 
iKobulnickv fc Dickevlll999l ) o r optical spectra tha t covered 
the Ca II absorption doublet (jSmoker et al.ll200"5l ). On the 
other hand, most important ions (e.g., Mg i, Mg ii. Si ii. 
Si IV, and C iv) have their resonance transitions in the 
rest-frame ultra-violet (UV) region. Their observed wave- 
lengths are still in the UV region for low-redshift targets 
like the MB. 

In this paper, we report on the general physical condi- 
tions and chemical composition of the MB gas far from 
the location of stellar populations, through quasar ab- 
sorption line analysis. We chose the sight-line toward a 
radio-loud quasar PKS 0312—770 at z = 0.223 that passes 
through the MB (Figure 1) at a heliocentric distance of 
~50 kpc (similar to the distance of the SMC). Because 
the quasar is radio loud we are also able to detect neu- 
tral h ydrogen gas as H i 21 - cm em ission/absorption fea- 
tures. IKobulnickv fc DickevI (|1999D reported detections of 
such features with three components at a heliocentric ve- 
locity of Vq ~200 km s^^. Based on their fitting results, 
they estimated a total neutral hydrogen column density 
of Nhi = 1.2x10^° cm~^ and mean spin temperatures 
of (T5) = 22, 29, and 46 K for the three components. 
Although the measured spin temperatures are contam- 
inated by warm gas structures that are not related to 
the absorbing components, they are simil ar to the spin 
temperatur es of the LM C (Tg = 30-40 K; iMebold et all 
1997: Dicke y et al.l ll994D and the SMC (Ts = 20-50 K; 
Dickey et al. 1999), and slightly small er than that of M31 
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Smoker et al.l (|2005[ ) detected Ca 11 K lines 



{Ts ^70 K; 
Recently, 

from the MB gas in their optical, medium-resolution spec- 
trum of PKS 0312—770. Here we report on detections of 
30 additional metal absorption lines in an HST/STIS UV 
high-resolution spectrum of the same quasar. We report 
the results of photoionization models in order to probe the 
physical conditions in the MB along this sightlinc. 

We describe the observations and data reduction in § 2 
and our procedures for line pr ofile fitting and ph otoion- 
ization modeling with Cloudy (jFerland et al.llT998l ) in § 3 
and § 4. We present our results in § 5, and discuss them 
in § 6. 

2. OBSERVATIONS AND DATA REDUCTION 

A UV spectrum of the quasar, PKS 0312-770 (my 
— 16.2), was taken with the Hubble Space Telescope 
(HST)/ Space Telescope Imaging Spectrograph (STIS) 

^ It is typically ~6 because only a tiny range of wavelength around Lya 

* IRAF is distributed by the National Optical Astronomy Observatories 
Science Foundation. 



(|Kobulnickvll2000[ ). using two echelle gratings: the E140M 
grating covering 1150A-1700A and the E230M grating set- 
ting that covers 2130A-2985A, with 6-10 km s'^ velocity 
resolution. The total exposure times are 37,908 sec for 
the E140M grating and 6060 sec for the E230M grating, 
which results in a spectrum with a S/N ratio of 6-28 per 
pixel] (AA ~0.015A) at A <1700A and 3-8 per pixel (AA 
'-~^0.04A) at A >2130A. A journal of the observations is 
shown in Table 1. 

The data were pro cessed with the sta ndard HST/STIS 
pipeline, CALSTIS ( Brown et ani2002f ). Continuum fits 
were made with standard techniques fChurchiU fc Vog^ 
I2OOID using the IRAF SFIT taskQ The normalized spec- 
trum, after rebinning every O.lsAfor A <1700Aand every 
O.40A for A >2130A, is presented in Figure 2 along with 
the 1(7 error spectrum. 

The STIS spectrum covers various absorption lines de- 
tected at > 5a, including (Lya, Mg i, C i, O i, N i, Mg 11, 
Mn II, Fe 11, Si 11, Ni 11, S 11, C 11, Si iii. Si iv, C iv, and 
N v) arising either from local absorbers in the Milky Way, 
from high- velocity clouds, from the MB or from higher red- 
shift absorbers a t z = 0.1983, 0.2018, 0.2026, and 0.2029 
(|Giandom 1 12005[ ) . Absorption lines arising from the MB 
(that are our main targets) are marked in Figure 2, and 
their rest-frame equivalent widths or detection limits are 
listed in Table 2. We can effectively separate absorption 
profiles of the MB from those of the Milky Way and the 
high- velocity clouds, because their velocity separations are 
large enough (the MB and the HVC are redshifted from 
the Milky Way by ~200 and 300 km s"\ respectively). 
Only Lya is not separated because it is hidden in the 
broad damping wing of the Milky Way absorption. The to- 
tal column density of this system (including contributions 
from the Milky Way, the MB, and the HVC altogether) 
was m easured to be log (H i/[cm~^ ]) ~20.85 (|Giandoni I 
[2005I) . or -20.87 (|Lehner et al.ll200l . 

As estimated below using photoionization modeling, the 
total H I column density of the MB toward the quasar 
PKS 0312-770 is logA^ (H i/[cm-^]) -20.2, just below 
the lower-limit of the criterion for damped Lya system 
(DLA). The total rest-frame equivalent width of Mg 11 
(VP;(2796) = I.91A) is also large enough for classifica- 
tion as a strong Mg 11 system (whose minimum equiv- 
alent width is Wr(2796) — 0.3A). Absorption systems, 
classified as strong Mg 11 systems, are usually associated 
with bright (L > 05L* ) galaxies within 40h~^ kpc (e.g., 
iBergeron fc Boiss^ll99lt) . The ionization state of the sys- 
tem, in which we detect very strong low-ionization tran- 
sitions such as O I and Fe 11 and which has either weak 
or undetected high- ionization transitions (e.g., C iv and 
N v), is very low, a s is often the case with DLA systems 
(e.g.. lLu et al.lll99^ . Thus, this system in the MB is prob- 
ably a local counterpart of DLA systems at higher redshift. 
We also confirmed that all absorption features in this sys- 
tem are clean without contamination by absorption lines 
of other systems at higher redshift. 

Visual inspection of the velocity plot of the MB sys- 
tem also suggests that the physical and chemical condi- 
tions of this system are somewhat similar to those of an- 

emission line has higher S/N ratio, 
which are operated by AURA, Inc., under contract to the National 
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other absorber detected i n the MB toward the young stars, 
pi 13 88 and DGIK 975 (iLehner et al.ll200lL 120081 : iLehned 
|2002|) : (i) having absorption lines with various ionization 
potentials from O i and C i to C iv and Si iv, and (ii) 
showing stro n g neu tral transitions such as O i and N i. 
iLehner et al.l ( 200l[ ) derived several properties of the sys- 
tem toward DI 1388: (a) dust depletion is moderate and its 
pattern is similar to that of the Galactic halo, and (b) the 
metallicity is lower than that in the solar neighborhood by 
1.1 dex, and still lower by 0.5 dex than in the SMC. Using 
photoionization models, we estimate these parameters to- 
ward the quasar, PKS 0312—770, separated by ~4.1 degree 
(i.e., ~4.3 kpc at the distance of the SMC) from DI 1388, 
so that we can study the variety of physical conditions in 
the MB. 

3. LINE PROFILE FITTING 

We determine the number of components needed to 
reproduce the observed spectrum and measure their 
line parameters (i.e., column density, Doppler param- 
eter, and radial velocity) using Voigt profile fitting. 
iKobulnickv k. DickevI (|1999f) detected three H i 21-cm ab- 
sorption components in their radio spectrum of PKS 0312- 
770 taken with the Australia Telescope Compact Array 
(ATCA). The H i 21-cm absorption lines are due to cool 
atomic hydrogen regions that probably have similar phys- 
ical conditio ns to DLA systems. 

Following IKobulnickv fc DickevI (|l999f ). we first used 
three H i components to fit other metal absorption pro- 
files in our HST /STIS spectrum. However, we found a 
three component fit insufficient bec ause the line widths of 
the th ree components presented in IKobulnickv fc DickevI 
(fT999l) (FWHM = 9.6, 13.2, and 10.5 km s-\ correspond- 
ing to &(H) = 5.8, 7.6, and 6.3 km s"-'^) could not reproduce 
the wing profiles at both sides of the metal absorption fea- 
tures. For example. Figure 3 shows a comparison between 
the observed and modeled O i A1302 profiles. We assumed 
that the Doppler parameters of oxygen and hydrogen were 
the same for this purpose, i.e. that turbulence dominates 
6, since that yields the maximum value for 6(0). If we 
increase the column density of O i, damping wings appear 
before our model reproduces both sides of the observed 
spectrum. To resolve this discrepancy, we refit the H i 21- 
cm absorpt ion lines, using the V oigt profile fitting code 
(MINFIT; IChurchill et al.ll2003l) (instead of using Gaus- 
sian fits as in [Kobulnicky & Dickey (1999)). 

These measurements of the H i column densities require 
some assumptions. We calculated them by 



sub-DLAs (i.e., logA^ (H i/[cm-2]) > 19)). The total 
H I column density after adding up the three compo- 
nents will be A'h I — 1.27x10^° [cm^^], which is consis- 
tent with the values measured in past studies using dif- 
ferent radio telesco pes with different beams on the sky: 
1-5x10^ ° fcm-^1 (iMathewson fc FordI I1984D. 1-2x102° 
cm-^l (IKobulnickv fc DickevI I1999D. 1.7 x10^° fcm-^l 
Kal berla et all I2005D . 1.3x10^° [cm-^] (|Lehner et al.1 
,20081 . and 7.1x10^° [cm"^] as an upper limit i ncluding 
contributions from our Galaxy (|Giandoni ll2005l) . These 
H I line wi dths are slightly broader tha n those found pre- 
viously by IKobulnickv fc DickevI (|l999t ) . Using three com- 
ponents with these larger Doppler parameters to fit the 
O I, we found it possible to reproduce the high-velocity 
side of O I A1302 by adjusting its column density to 
log (O i/[cm~2]) ~16, as shown in Figure 43- 

The only remaining disagreement is an underproduction 
by the model of O i A1302 at the low-velocity side of the 
absorption profile. Because a similar discrepancy is seen 
in other metal absorption lines, such as Mg ii, Fe ii, and 
Si II, there is probably an additional component at the 
low-velocity side of this system. Assuming a fiducial line 
width of & = 10 km s~^ (i-C, the average Doppler param- 
eter of the other three H i 21-cm lines is 9.0 km s~^), we 
adjust both the line position and column density of O i 
and find the best combination: Vq = 160.8 km s^^ and 
logA^ (O i/[cm-2]) ^15. 

In fact, another strong constraint for this additional 
component comes from the radio observation. We can 
place an upper limit on the H i 21-cm line optical depth 
of this additional component at the line center, because 
we did not detect it at the appropriate wavelength in the 
H I 21-cm absorption spectrum. The S/N-ratio of the 
spectrum is ^40 at that wavelength, which corresponds to 
a la optical depth limit of Tmax 



ln(l - a) -0.02 
As described in the next section, our best photoionization 
model for this component gives log (H i/[cm^2]) of 19.31 
and 6 of 13.0 km s^^ for the H i line, whose optical depth 
at line center would be r = 0.022, smaller than the Icr de- 
tection limit in the radio spectrum (Figure 5). We assume 
the same spin temperature for this component as was mea- 
sured for the closest component at Vq = 175.8 km s~^ (Tg 
~ 22 K). The line parameters for this additional compo- 
nent, determined from a fit to the O i are listed in the 
fourth line of Table 3. However, we will not consider this 
component in the discussion (§ 6) because the fitted line 
parameters are all based on an arbitrary Doppler param- 
eter. 

4. PHOTOIONIZATION MODELING 

4.1. Modeling Procedure 

We briefly summarize our modeling procedure, although 
it is similar to previous studies (e.g., Churchill & Charl- 
ton 1999) . Using the photoion ization code Cloudy, version 
07.02.00 (jFerland et al.lll998D . for each of the five clouds 
(the fifth one is introduced later) we search for the best 
combinations of the fit parameters: 1) ionization param- 
eter (log [/ — log[n^/ri/f], defined as the ratio of ionizing 

^ We also cannot reject the alternative solution that instead of a broader H I component on the high-velocity side of the system, there are two 
components blended together at that velocity. Such a model, however, would involve three additional free parameters (i.e., position, column 
density, and line width) for which we do not have constraints from our radio spectrum. We adopt the one component model for simplicity. 

i.e., a typical uncertainty of our estimates of H I column densities is also about 2.5%. 



Nhi = 1.823 X 10 
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/ 



T2ldV, 



(1) 



where Ts is the spin temperature in degrees Kelvin, 
and r2i is the H i 21-cm optical depth. We adopte d 
the spin temperatures from ' Kobulnickv &: DickevI (jl999l l. 
Then by integrating the optical depth of the model pro- 
file that we best fit to the observed spectrum, we cal- 
culated the H I column densities. Our best fit param- 
eters are listed in the first three rows of Table 3. The 
column densities are large enough to be classified as 
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photons to the number density of hydrogen in the absorb- 
ing gas) and 2) metaUicity (log(2'/Zo)) in solar unitsF^ 
The MB is hkely to be confined and compressed by its in- 
teraction with th e Galactic halo, but this does not produce 
shock ionization (|Bland-Hawthorn et al.|[2"007( ). Therefore, 
we consider only photoionization. We optimize on the ob- 
served column densities of H i or other metal ions (i.e., O i 
and Si iv) that are listed in the first column of Table 3, that 
is we require the Cloudy models to produce those column 
densities. For example, in the case of the component at 
Vq ~176 km s~^, we always fix the H i column density to 
log TV (H i/[cm~^]) = 19.59, and search for the best values 
of logJJ and \og{Z/ZQ)m a grid with steps of 0.1 dex, to 
reproduce as many absorption lines from other transitions 
as possible. We repeat this procedure for the other four 
components, individually. We assume that the absorbers 
arc plane-parallel structures of constant density, and that 
they are in photoionization equilibrium. We initially adopt 
a solar abundance pattern, but also explore some varia- 
tions based on the observational constraints. For the in- 
cident radiation field, we first consider a pure extragalac- 
tic background radiation with contributions from quasars 
and star forming galaxies (with a photon escap e fraction 
of 0.1), following iHaardt fc Madaul (1 19961 . [200111 . We also 
explore the effects of other incident radiation fields, includ- 
ing fiuxes from our Galaxy and the LMC, as described in 
the next section. 

For a given logJ7 and log(Z/ZQ), Cloudy computes 
the column density of each element in various ionization 
stages, and the equilibrium gas temperature, T. We can 
calculate the Doppler parameter for each element using b 
~ + ^Lrb' where br is the thermal broadening de- 
fined as 6t = y^2kT/m and bturb is the broadening from 
gas turbulence and bulk motion. We use the measured 
b for the transition on which we optimized (listed in Ta- 
ble 3) in order to calculate bturb, which is then applied for 
other elements. Using these derived line parameters from 
the model, we synthesize a spectrum, after convolving the 
instrumental line spread function of HST/STIS, and com- 
pare it to the observed spectrum. The synthesized spectra 
are compared to the observed spectrum by eye, because 
models with metallicity and ionization parameters that 
differ only slightly from the best values (e.g. 0.1 or 0.2 
dex) would deviate si gnificantly from the observed spec- 
trum (see Figure 4 of iMisawa et al.l [20081 ) . Moreover, a 
formal procedure like fitting is not applicable because 
the system in the MB has various (more than 10) transi- 
tions that we must consider simultaneously, and it is very 
difficult to deter mine how they shoul d be weighted in a 
calculation fe.g.. iMisawa et al]l2008f) . 

The transitions we optimize are H i and O i, and both 
have very low ionization potentials (IP = 13.6 eV). As 
frequently reported for photoionization models of Mg ii 
absorbers, the low-ionization phase clouds that produce 
Mg II absorption lines also produce other low-ionization 
transitions (e.g., C ii. Si ii, and Fe ii), but not high- 
ionization transitions (e.g.. Si iv and C iv). An addi- 



tional high-ionization phase is almost always required to 
reproduce these transitions. Therefore, we repeat line fit- 
ting and photoionization modeling for the high-ionization 
phase by optimizing high-ionization transitions (Si iv in 
this study, as described below). Finally, we synthesize a 
model spectrum including both low- and high-ionization 
phases, and compare it to the observed spectrum. 

As described above, our method for deriving constraints 
on metallicities and ionization parameters of the absorbing 
gas relies on Voigt profile fitting of separate components 
of the absorption profile of a particular transition which 
we feel is best constrained. We then use photoionization 
modeling to infer the column densities and Doppler param- 
eters of other transitions in these Voigt profile component 
clouds. Synthesized profiles from these models are com- 
pared to the data in order that we can place constraints 
on the parameters. Although this method does rely on 
the assumptions behind the Cloudy models, it has some 
advantages over simply taking the ratios of apparent col- 
umn densities of selected transitions in order to determine 
metallicity. First of all, we can separately examine prop- 
erties of individual clouds along the line of sight at dif- 
ferent velocities. We are not averaging these components 
together, which the apparent column density method re- 
quires since it does not distinguish, for example, how much 
of the hydrogen is associated with the separate clouds. 
Secondly, we can use the components determined from un- 
saturated lines in order to constrain the properties of those 
clouds using other saturated or partially saturated compo- 
nents as well. The measured Doppler parameters for the 
unsaturated lines can be used along with the temperature 
given by a given Cloudy model in order to determine model 
Doppler parameters of other lines that may be saturated 
or blended. Comparing the model to the shapes of the ob- 
served profiles of these lines, particularly the shapes of the 
sides of the profiles, often yields meaningful constraints on 
parameters. Models also take into account the appropri- 
ate ionization corrections in each case. In this way we can 
determine the range of acceptable parameters, consider- 
ing all observational constraints. Finally, we are able to 
consider separate phases of gas, having different densities 
and velocities along the line of sight, by comparing model 
predictions to the observations, component by component. 

4.2. Alternate Incident Radiation Fields 

Because the Milky Way (MW) and the LMC are located 
within several tens of kilo-parsecs of the MB, they also con- 
tribute as additional ionizing radiation sources. Therefore, 
we consider three alternative inci dent radiation fields, (i 
extra-galactic background fEGB: [Haardt fc Madaul [T99i 
l2001h plus radiation fr om the MW at a distance of -D = 
50 kpc from the MW (iFox et al.ll2005a[) . (ii) a maximum 
flux model, with the EGB, the MW radiation, plus radia- 
tion from the LMC with a 30 % escape fraction, and (iii) an 
intermediate case, with EGB, the MW radiation, plus ra- 
diation from the LMC with a 15 % escape fraction. We 
constr uct thes e radiation fiel ds following the procedure of 
iBland-Hawthorn fc Malonevi (1999 , 2002) . Figure 6 shows 

The default solar composition we assume is listed in Hazy 1, a manual of Cloudy. For example, [Si/H]0, [N/HJq, and [0/H]q are —4.46, 
—4.07, and —3.31, respectively. 

We do not consider the contribution from local O/B-type stars in the MB, because it is less likely that our sightline to the background 
quasar goes through stellar associations, at least c ompared to the sightlincs toward stars in the MB. It is also suggested that O/B-type stars 
tend to localize in the wing of the SMC l llrwin et al.iil99a: .Battinelli Demers„1992iV 
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the strength of the ionizing radiation from the MW and 
the LMC as a function of a distance from the center of 
the MW. At D -50 kpc, at which the MB is located, the 
contribution from the LMC to the radiation field is signif- 
icant. On the other hand, the contribution from the SMC 
is negligible because of its low luminosity. 

The spectral shapes for each of the four incident radia- 
tion field models on the MB gas are plotted in Figure 7. 
The radiation from the MW and the LMC start to dom- 
inate the EBR at log (i^/ [Hz]) < 16.1, and strongly domi- 
nate over the EBR at log(zy/[Hz]) < 15.5. Normalizations 
for these fields, in units of the density of ionizing photons, 
and the contribution from the local flux sources (MW and 
LMC) compared to the EBR are summarized in Table 4. 

5. RESULTS 

Transitions that are detected in the HST/STIS spec- 
trum and those that provide useful limits are shown in 
Figure 8, from the lowest (Mg i, IP = 7.6 eV) to high- 
est (N V, IP = 97.9 eV) ionization potentials, following 
the Lya profile in the first panel. We also display the ob- 
serve d H I 21-cm absorption profile (|Kobulnickv fc DickevI 
[T99I . InFi gure 8, km s ^ denotes the apparent optical 
depth-weighted median of the MB absorption system, cor- 
responding to a heliocentric velocity of V© = 209.9 km s~^. 
A number of absorption lines from various ions (i.e., Mg i, 
C I, O I, N I, Mg II, Mn II, Fe 11, Si 11, Ni 11, S 11, Si iv, and 
C iv) are detected. For N v, we can place only an upper 
limit on the equivalent width. The Mg 11 A2803 profile 
suffers from a data defect on its blue side, so we use the 
Mg II A2796 profile as the main Mg 11 constraint. Since 
the Mg I A2853 profile is noisy (i.e. S/N ~3.7 per pixel), 
it is only used as a loose constraint on models. 

As a starting point for photoionization modeling, we 
roughly estimate possible ranges of gas temperature and 
ionization parameter (logt/). The observed line widths of 
all transitions are similar, which suggests that bulk motion 
(gas turbulence) is a dominant source of line broadening. 
The relatively narrow 21- cm lines, especially the one with 
b = 6.4 km s~^, then imply that the gas temperature is 
very low. Such a low temperature implies that the ion- 
ization parameter is also very small {\ogU < —5.0). For 
the metallicity, the value in the SMC has been estimate d 
asZ = 0.1 - 0.2 Zq (iPagel et al.lll978t IWeltv et al.ll2001h . 
Therefore, in the following we explore metallicities and 
ionization parameters, in steps of 0.1 dex, in the ranges 
log {7 = -7.0 - -4.0 and log(Z/Zo) = -1.0 - 0.0 unless 
other parameter ranges are suggested. 

5.1. Three Hi 21-cm Clouds 

At first we seek the best photoionization model param- 
eters for the three H i 21-cm clouds at Av = —34, —3, and 
26 km s~^ from the system center (i.e., V© ~ 176, 207, and 
236 kms^i) that we determined in § 3, by optimizing a fit 
to the H I 21-cm absorption line. We require our model to 
reproduce the column density of H i from that fit. Because 
the cloud at Av — —3 km s~^ has clear detections in C i, 
N I, Ni II, and S 11, without blending with other absorp- 
tion features, we begin with this cloud. We can constrain 
the metallicity to be -0.9 < \og{Z/ZQ) < -0.5/-0.5 < 
\og{Z/ZQ) < —0.2 to avoid over/under-production of Ni 11 
and S II. Therefore, the acceptable metallicity is \og{Z/ZQ) 



'^— 0.5. With this metallicity, the ionization parameter 
must be —6.0 < \ogU < —4.8 to avoid over/under pro- 
duction of C I. Even our favored model, with log(Z/Zo) 
'-~^— 0.5 and logC/ ~— 5.8, for this Aw = —3 km s~^ cloud 
substantially overproduces the N i A1201 absorption. To 
reconcile the model with the observed N i A1201 absorp- 
tion, we decrease the nitrogen abundance by a factor of 0.7 
dex compared to its solar abundance pattern. The pro- 
duction mechanisms of nitrogen are poorly understood. 
However, a nitroge n deficiency in the S MC has already 
been reported (e.g.. lMallouris et"ani200l[ ). a nd similar de- 
ficien cies are seen in some DLA systems (e.g., Pettini et al.l 
I2002D and weak Mg 11 systems (|Zonak et al.l l20 04 !) . Since 
O I and H i are strongly coupled in terms of their ion- 
izatio n, their ratio is oft en used as a measure of metal- 
licity (|Lehner et al.l [20081 ). In our case, the O i is highly 
saturated so that the Voigt Profile components cannot be 
measured directly from the profile. However, once we have 
determined the metallicity from the weaker Ni 11 and S 11 
profiles, we can verify that the favored model produces 
an O i/H I ratio consistent with these metallicities. For 
the —3 km s~^ component we find log[N(0 i)/N(H i)] = 
—3.81, corresponding to log{Z/ZQ)= —0.5. The fact that 
this is consistent with the metallicity inferred from the 
weaker profiles confirms that ionization corrections are in- 
ferred properly from our models, and that the saturated 
O I profile has been separated self-consistently into Voigt 
profile components. 

Next, we consider the H i cloud at Av — 26 km s^^. Be- 
cause this component has the largest recessional velocity, 
it should reproduce at the high- velocity side of all absorp- 
tion features seen in the O i, Mg 11, Fe 11, and Si 11 lines. 
We also require that the component reproduces the weak 
Ni II absorption. The Si 11 provides the best lower limit on 
metallicity, \og{Z/ZQ) > —0.8 because the less saturated 
Si II A 1304 profile is available. Similarly, an upper limit 
of log(Z/Zo)< —0.6 applies in order that Ni 11 is not over- 
produced. At the preferred value of log(Z/Zo) ~— 0.7, 
there is a strict lower limit on the ionization parameter of 
log [/ > —6.1 under which C i would be over-produced. 
We do not place any formal upper limit on the ionization 
parameter of this cloud, since the constraints from the rel- 
atively low signal-to-noise spectrum were not significant. 
However, we can place a marginal upper limit, logU < 
—5.0 to avoid an under-production of Si 11, once we adopt 
\og{Z/ZQ) — —0.7. We favor values toward the lower end 
of this range, \ogU ^—6.0, because the fits to O i, Mg 11, 
and Fe 11 are slightly better. We also find that this cloud 
must have a deficiency of nitrogen of 1.4 dex compared to 
the solar value, although this value could be smaller (§. 6). 

Because the third H i cloud at Av ~ —34 km s^^ is 
heavily blended with the adjacent clouds, it is difficult to 
place strict constraints on its physical conditions. There- 
fore, we simply note that its metallicity and ionization 
parameter could be very close to those of the cloud at Av 
= 26 km s^^. Several transitions, Mg i, Ni 11, and C i, are 
over-produced for metallicities higher than \og{Z/ZQ) > 
—0.6 or —0.7. The metallicity should also be greater than 
— 1.0, to avoid an under-production of Si 11. Although we 
cannot place strong constraints on the ionization param- 
eter, the \ogU should be between —6.0 and —5.0 in the 
acceptable range of the metalHcity, — 1.0 < \og{Z/ZQ) < 
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—0.7, to reproduce the observed profiles of C i, Si ii, and 
O I. We adopt a model with logU '^—5.7 and \og{Z/ZQ) 
'~— 0.7, the same as for the Av = 26 km s~^ cloud, because 
it is consistent with the observations. We also find that 
this cloud must have a deficiency of nitrogen of 1.0 dex 
compared to the solar value. 

5.2. High Ionization Cloud 

In addition to the three H i clouds above, one high- 
ionization cloud, at Av = 7kms~^ (i.e., Vq = 217km s~^), 
is necessary to reproduce absorption in the high-ionization 
lines. Si iv and C iv. To determine the photoionization 
model parameters for this cloud, we optimizes on the Si iv 
column density, since the Voigt profile fitting is better in 
this region than it is for the C iv doublet. The result of 
a Voigt profile fit to the Si iv is given in the last row of 
Table 3. Comparing to the observed C iv absorption, we 
find an ionization parameter of —2.7 < log?7 < —2.4. The 
metallicity of the high-ionization cloud is constrained to 
be log(Z/ZQ) > —4.0 by the requirement that the corre- 
sponding H I 21-cm line is not detected at that velocity. 
The metallicity could actually be much higher, as high as 
that of the low ionization clouds, but we have no way to 
place further constraints. 

5.3. Additional i Cloud 

A model with three low ionization clouds and one high- 
ionization cloud reproduces the observed spectrum very 
well except for the low-velocity side of O i, Fe ii, and 
Si II. As noted in § 3, we add an additional cloud at Av 
= —49 km s~^ (i.e., Vq = 161 km s~^) in order to fully 
produce the observed O i absorption. For this cloud, we 
can constrain the metallicity by the requirement that the 
optical depth of H i 21-cm at line center must not ex- 
ceed T 0.025. We obtain \og{Z/ZQ) > -1.0 for logC/ > 
-5.1, and \og{Z/ZQ) > -0.9 for logC/ < -5.2. The ion- 
ization parameter should be logf/ > —6.0, below which 
Mg I would be overproduced. An upper limit on the ion- 
ization parameter of this component is —3.5, to avoid an 
over-production of Si IV. We again see that nitrogen is de- 
ficient in this cloud by 0.6 dex. We give a sample of an 
acceptable model for this cloud, with \og{Z/ZQ) = — 1.0 
and logC/ = —5.1 in Table 3. However, we note that the 
parameters could instead be more similar to those of the 
other clouds. 

5.4. Our Best Model 

Table 3 lists the best parameters of our photoionization 
models: column (1) is the optimized transition, column (2) 
is the heliocentric velocity of the absorption line, column 

(3) is relative velocity from the system center, columns 

(4) and (5) are absorption optical depth at the line cen- 
ter and spin temperature measured from the H i 21-cm 
emission line (only for H i 21-cm lines), columns (6) and 
(7) are the best Voigt profile fit values of column density 
and Doppler parameter, columns (8) and (9) are the best 
model parameters (assuming a solar abundance pattern) 
of metallicity and ionization parameter, column (10) is gas 
temperature from our model, columns (11) and (12) are 
the gas volume density and the thickness of the absorber 
assuming a plane-parallel structure, and column (13) notes 



changes in the nitrogen abundance (compared to the So- 
lar abundance pattern) that arc necessary to reproduce 
the observed spectrum. The synthesized spectrum, using 
the best model parameters, is overplotted on the observed 
spectrum for each transition in Figure 8. If we compare the 
models with/without nitrogen deficiency to the observed 
spectrum, it is obvious that the nitrogen deficiency is nec- 
essary, as shown in Figure 9. A summary of how specific 
transitions were used to constrain log U and \og{Z/ZQ) for 
the MB system is given in Table 5. We also compare the 
equivalent widths of the observed and the modeled spectra 
in Table 2, and confirmed they are in good agreement. 

5.5. Effects from. Additional Radiation Sources 

Because this absorber is located in the MB, close to 
the MW and the LMC, the radiation field around it could 
be enhanced by contributions from these galaxies, as men- 
tioned in § 4.2. We perform photoionization analysis using 
the various radiation fields listed in Table 4. Because the 
purpose of this analysis is to examine how the shape of ra- 
diation field affects the results of the photoionization mod- 
eling, we always use the best- fit parameters from Table 3. 
We compare the model spectra using various incident radi- 
ations to the observed spectrum in Figure 10. On the ob- 
served spectrum, we do not see any significant differences 
between the results for the different radiation field models, 
except for C i A1329, for which the model is slightly im- 
proved if we take a radiation from the LMC into account. 
However, the quality of spectrum around this line is not 
very high. Apparently the same fit parameters are still ac- 
ceptable for the alternative radiation fields. However, once 
we estimate the line parameters of various transitions in 
the strongest component at Av = —3 km s^^ when us- 
ing these incident fields, there are noticeable differences 
for some transitions (i.e., C iv, N iii, Al iii, and Si iv) 
that are very weak and/or positioned at low-S/N regions. 
This is because the number of hydrogen ionizing photon is 
basically same by our assumption (because we assume the 
same ionization parameter), however, the ionizing photons 
of some ions (especially ions with low ionization potentials) 
would be increased significantly once the additional radia- 
tion is considered. Although the ionization parameter that 
we infer is the same for these alternative radiation fields, 
the number of ionizing photons is substantially increased 
so that this same ionization parameter corresponds to a 
higher density. For the Av = —3 km s^^ component the 
density has increased from 0.56 cm~'^ to 19 cm~^, and the 
size decreased from 25 pc to 0.72 pc in the most extreme 
case. The other clouds arc similarly smaller, with the Si iv 
cloud reduced in size to 0.57 kpc under the more intense 
LMC radiation field. Since the differences in the ionization 
parameters that we infer are so small between the different 
models, it is not practical to distinguish between the dif- 
ferent possible radiation fields using these data, and thus 
the inferred densities and sizes are uncertain. Hereafter, 
we discuss the results using only the EGB radiation, as 
summarized in Table 3. 

6. DISCUSSION 

6.1. Origin of Metal-enrichment in the MB 

There have been three scenarios prosed to explain the 
current chemical composition of the MB. Tidal stripping 
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from the SMC, purported t o explain the origin of the MB 
([Gardiner &: Noguchilll996l ) would imply a chemical com- 
position similar to the SMC itself, or perhaps slightly more 
metal poor if the MB is preferentially formed from mate- 
rial drawn from the comparatively pristine outskirts of the 
SMC. The presence of any chemical gradient in the SMC is 
not ye t well defined, but the recent results of lCarrera et al.l 
(|2008f ) suggest it may be non-negligible. A galactic wind 
from the SMC could enrich the MB with metals, and 
such a wind would likely be enhanced in alpha elements, 
given the composition of supernova-driven winds in other 
dwarf galaxies where such signatures are observed several 
kpc f rom their origin (Marti n et al.ll2002t [Stricklan d et al] 
120041) . The final poss ibility is in-situ enrichment from 
stars formed in the MB (lOemers fc Battinellilll998D. Give n 
the ~200 Myr age of the MB (|Gardiner fc Noguchill 19961) . 
there has been ample time for several generations of the 
most massive B stars to evolve and contribute their nucle- 
osynthetic products to their surroundings. Some combi- 
nation of these processes are likely to play a role, but the 
b alance of these h as not y et been determined. 

iRolleston et al.l (119991 ) found an underabundance 
(~0.5 dex) in the light metals of a B-type star in the MB, 
DGIK 975, compared to the SMC H ii regions. Similarly, 
metallicity of the interstellar medium in the MB toward a 
young star, DI 1388, was measured to be ^0 .2 dex lower 
than the SMC H ii regions (jLehner et al.ll2008D . The sight- 
line toward our target, PKS 0312-770, has an angular sepa- 
ration of ~4.1° from DI 1388, which corresponds to about 
~4.4 kpc in physical scale at the distance of the SMC 
{d ~60 kpc). However, we have found that the metallic- 
ity in the MB toward PKS 0312-770 is higher than that 
measured toward DI 1388 by ~0.5 dex, and even ~0.1 
- 0.2 dex more metal-rich than the SMC H ii regions. 
What is the source of this difference? One possible idea is 
that the sightline toward PKS 0312-770 is not mixed with 
metal-poor gas, like the sightline toward DI 1388. If this 
is the case, absorbers with lower-metallicities should have 
higher total hydrogen column densities. Such a trend (i.e., 
a gradual increase in metallicity with decreasing hydrogen 
column density) has already been pointed out for various 
objects at higher redshift, including galaxies , (sub-)DLA 



systems, and the intergalactic medium (e.g . . iBoisse et al. 



2008). 
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There are observational results that support the sce- 
nario above. Column densities of sulfur (a modestly dust- 
depleted element) are similar between these two sight- 
Hnes. Each of the three clouds toward PKS 0312-770 (with 
logTV (S/[cm-2])is 14.2, 14.4, and 14.20 has a similar 
column density to the total (logiV (S/[cm^^]) = 14.35) 
in DI 1388. However, the PKS 0312-770 clouds have to- 
tal hydrogen column densities (logiV (H/[cm~^]) = 19.6, 
19.6, 19.70 0.6 dex smaller than the total in DI 1388 
(logA^ (H/[cm~^]) = 20.2). This imphes a higher metal- 
licity and less dilution in PKS 0312-770. These results are 



consistent with the scenario above: i.e., metal-enriched gas 
flowed out toward the MB from the SMC almost isotrop- 
ically, but the MB absorbers toward PKS 0312-770 were 
only weakly diluted by metal-free material, while the ab- 
sorber toward DI 1388 is significantly diluted. 

6.2. Nitrogen Deficiency in the MB 

The overall metallicity and the chemical abundance pat- 
tern in the MB can be used to infer th e ori gin of this inter- 
Cloud material. iHamblv et al.l (jl994j) and IRolleston et"al] 
(|l999) measured He, O, Si, Mg, and N abundances for 
several early B stars in the MB, assuming that their com- 
position reflects the present-day ISM, and concluded that 
these stars were 0.5 dex more metal-poor than SMC stars. 
These authors suggested that the MB stars formed from 
a mixture of SMC and un-enriched gas. iDufouil (| 19841 ) 
analyzed the H ii region He, C, N, O, and Si abun- 
dances for the SMC, finding that the overall metallicity, 
as measured by the alpha elements O and Si, is ~0.5 
dex more metal poor than the Sun. The N/Si ratio in 
H II regions is 0.5 dex lower than the Sun, consistent with 
the t rend observed in other low-metallicity dwarf ga laxies 
fe.g.. lKobulnickv fc Skillnm]ll996HNava et al.ll2006D . We 
summarize these measurements from the literature along 
with ours in Table 6. Locations of the MB stars are also 
shown in Figure 1. 

Figure 11 shows the abundance ratio log(N/Si) versus 
log(Si/H) for the aforementioned measurements from Ta- 
ble 6, assuming a solar abundance pattern if another ele- 
ment (instead of silicon) is used to estimate these param- 
eters. The figure includes our three velo city components 
towar d PKS 0312-770 (open circles), the IRolleston et al.l 
(jl999f) M B stars (filled s tars), the gas toward the MB star 
DI 1388 (jLehner et al.l (|200 8): open star), the SMC star 
AV 304 (open square), the Sun (solar symbol), the range 
of SMC H I I regio ns (filled square), the dwarf galaxies from 
iNava et al] (|2006f ) (crosses; values computed from O abun- 
dance measurements assuming a sol ar Si/0 ratio), and 
damp ed Lyman alpha systems from iHenrv fc Prochaskal 
(|2007l ) (dots). The abscissa reflects the overall alpha el- 
ement abundance of the systems in question and shows 
that the metaUicities of the SMC star AV 304, the MB 
stars, and the MB gas toward PKS 0312-770 are gener- 
ally consistent with the range of metaUicities in SMC H ii 
regions. The MB star DGIK 975 is a possible exception, 
appearing ^0.5 dex more m etal poor than th e rest of these 
measurements. Figure 3 of IRolleston et al.l (|1999D shows 
that this star also lies the furthest from the SMC, in a low- 
H I column density region equidistant between the SMC 
and LMC. 

The N/Si ratio shown on the ordinate of Figure 11 is a 
measure of the chemical enrichment timescale. Si is syn- 
thesized in massive stars and returned to the interstellar 
medium through supernovae on timescales of ~10 MvrF^ 
Nitrogen, by contrast, is thought to be produced in most 
galaxies by low- and intermediate-mass stars and released 



We fitted the H I 21-cm absorptio n profile with three components. The total S II column density (log(S II) ~14.75) is slightly lower than 
that mea sured inlLeh ner et ahl II2008I ) because we assume S II A1254 is blended with an unrelated line to avoid over-absorption of S II A1251, 
for which iLehner er^Ln2008l) applied the AOD method directly. 

Note that neutral hydrogen column densities toward DI 1388 and each of the three clouds toward PKS 0312—770 arc almost the same 
(log Nil I ~19-6), but the former is considerably more ionized (log Nn n ~20.0) than the latter (log A^h ii ~17.7). 

However, there is the possibility that a-elements release d by supernovae remain for an extended period in hot 10® K bubbles and require 
several hundred Myr to mix with ambient galactic material l|Tenorio-Taglail996i : iKobulnickv Skillmaniil997i) . 



8 



Misawa et al. 



on timescales of >100 Mvr (Ivan den Hoek fc GroenewegenI 
Il997t iMarieol [Mil iRenzini fc Volil ll98lM '"l Thus, the 
N/Si ratio may drop during prolonged starbursts and rise 
during prolonged periods of quiescence, functioning as a 
kind of "clock", marking the tir ne since the most recent 
major epi s ode o f star form ation (lEdmun ds fc Pagel 1978t 
Pantelakil 119881: iGarnettI [19 90: Kobuln ickv fc Skillmaij 



1998( ). iHenrv fc Prochaskal ~(2007) used this scenario to 



model the chemical evolution of DLA systems in the N/Si 
vs. Si/H plane and found that they could reproduce the 
properties of most DLA systems as a function of two vari- 
ables: the star formation efficiency and the age of the 
system. They concluded that small ages (i.e., less than 
250 Myr) are required to produce the low N/Si systems, 
while high star formation efficiencies lead to higher Si/H 
and lower N/Si ratios. However, a number of other pro- 
cesses may contribute to shaping the chemical evolution of 
a system. The arrows in Figure 11 show qualitatively the 
four possible evolutionary vectors attributed to N enrich- 
ment (vertically upward) , a enrichment (toward the lower 
right), dust depletion (toward upper left), and dilution by 
metal-poor gas (leftwards). 

Figure 11 shows that the 207 km s~^ component to- 
ward PKS 0312-770 is 0.3-0.4 dex more metal-rich than 
the bulk of the SMC and other MB stars. This compo- 
nent also exhibits an N/Si ratio that is 0.2 - 0.3 dex lower 
than the SMC stars and H ii regions. This departure is 
consistent with a composition consisting of SMC material 
augmented by a small amount of a enrichment from su- 
pernovae. None of the MB stars in Figure 11 nor the ISM 
probed toward the MB star DI 1388 share this abundance 
pattern. EH The PKS 0312-770 line of sight is therefore 
the first MB location with a metallicity that is similar to, 
or perhaps slightly higher than, the SMC at large. 

Figure 11 shows that the absorbing component at 
236 km s^^ has a metallicity similar to the SMC and most 
MB stars, but the N/Si ratio is about 1.0 dex lower [^ 
Such a low N/Si ratio places this material among the most 
N-deficient damped Lyman alpha systems and suggests a 
nucleosynthetic history dominated by a-pr oducing mas- 
sive stars. In the lHenrv fc Prochaskal (|2007D models, such 
a low N abundance at such high metallicity can only be 
achieved by very efficient and very recent star formation 
such that massive stars dominate the mass-averaged nu- 
cleosynthetic contribution with virtually no contribution 
from longer-lived N-produci ng stars. Given the p resence 
of B stars in the MB (e.g., iRoUeston et al.l [19991 ). some 
supernova activity during the < 200 Myr lifetime of the 
MB is likely. In situ enrichment appears a plausible expla- 
nation for this component on the basis of chemical cloud. 
However, it is also possible that a Si enrichment of ambient 
SMC material plus dilution from putative metal-poor halo 
gas could also explain this data point (i.e., a combination 
of vectors that together drive evolution in Figure 11 in a 
downward vertical direction from the SMC H ii regions 
composition) . 

Finally, the 176 km s~^ component has a Si/H ratio 



consistent with the SMC and the 236 km s~^ component, 
but the uncertainties on the N/Si ratio are too large to 
place meaningful constraints on the origin of this mate- 
rial. Moreover, abundance pattern would also be strongly 
affected by dust depletion and ionization conditions, which 
can only be explored with very high-S/N spectrum. There- 
fore, this component may either be N-deficient or it may 
be consistent with the SMC H ii regions. 

6.3. Dust Depletion in the MB 

In addition to various abs orption lines w e have detected 
in the HST/STIS spectrum, ISmoker et all (|2005t) also de- 
tected the Ca II K line in their medium-resolution (i? = 
6000) optic al spectra of 7 quasa rs behind the MB including 
our target ([Smoker et al.ll2005[ ). By comparing the total 
column densities of Ca ii and H i (measured from H i 21- 
cm emission line), they found that the abundance ratio 
of Ca II to H I in the MB is systematicall y higher than 
that o f Galactic gas by a factor of ~0.5 dex. I Smoker et all 
(|2005f ) proposed possible scenarios for this difference, such 
as the higher ionization condition of the hydrogen gas, and 
weaker dust-depletion in the MB. To test the scenario, a 
higher-resolution spectrum will be necessary to deblend an 
unresolved Ca ii K profile into multiple components (as we 
did for other transitions in the UV spectrum) in order to 
constrain the photoionization model. 

In § 5, we found the best model parameters. However, 
there is still an ambiguity due to the possible effects of 
dust, which would lead to different inferred physical con- 
ditions. Our data was not of sufficient quality to measure 
dust depletion, but the estimated low gas temperature of 
the three H i clouds {Tgas < 1000 K) would imply that 
dust grains could survive in the absorber. Moreover, the 
Ca II K absorption strength, expected from our best model 
{Wrest=0 -99 A), is seven times greater than the observed 
value toward PKS 032-770 fW;..,,=0.14A: ISmoker eHI] 
l2005f ). This also implies that the absorber contains sub- 
stantial amounts of dust be cause calcium is one of the most 
severely deple ted elem ents (Sav age fc Sembach[| 19961 ) . Ac- 
tually, Lehner et al.l ((2008 1 suggested Si and Fe in the 
MB toward DI 1388 are depleted to dust by factors of 
—0.45 and —0.61 dex respectively, although the condition 
could be different toward PKS 0312—770. If depletion onto 
dust is significant, then the Si abundances and the implied 
metallicities become even larger, making this sightline sig- 
nificantly more metal-enriched than the SMC itself. 

6.4. Comparison with Damped Lya Systems at high-z 

Damped Lya systems (DLAs) are characterized by high 
H I column densities (i.e., logA^ (H i/[cm~^]) > 2 x 
10^0 cm-3) and low metallicities (i.e., Z r^Q.l - 0.01 Zq) 
(e.g. iPettini et ai]|1997l : lProchaska et al.ll2003D . DLAs are 
also known to have lower nitrogen abundance relative to 
g-elements, compa red to solar abundance pattern (e.g., 
iPettini et al.l 12003 ). DLAs provide plentiful information 
on the physical condition in gas clouds in the ancient uni- 
verse that can never be traced by stellar objects. However, 



See, however. lAckerman et aTl II2004I ) : [Spite et al.l (120051 ) : [Pettini et al.l I I2008I) for evidence that massive stars may dominate N production 
at very low metallicities, Z < 0.01 Zq. 

Hot stars in the MB could be affected by additional production of nitrogen. 

Normalization of the spectrum may be performed incorrectly to the low-S/N (~5) spectrum region around N I lines. If we reduce nitrogen 
abundance by only —1.0 dex from the solar abundance pattern (i.e., same as the 176 km s-^ component), N I A1200b and N I A1201 would be 
reproduced better (see Figure 12). 
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so far DLA absorbers have not been identified clearly, espe- 
cially at higher redshift, alt hough at least ellipt ical galax- 
ies are probably ruled out (^ Calura et al.l [20031 ). Thus, it 
is quite helpful to find local counterparts of those systems 
and study their properties in detail. 

With respect to H i column density, the MB absorbers 
toward PKS 03f 2— 770, whose total column density is just 
below the criterion to be classified as a DLA system, could 
be analogs of high-z DLA systems. However, their nitro- 
gen abundance relative to a-elements (e.g., [N/Si]) tend 
to be small and positioned at the lowest end of the [N/Si] 
distribution of DLA systems in Figure 11. Particularly, 
the nitrogen abundance of the Vq = 236 cloud is very 
small compared to that expected for DLA systems and 
blue compact galaxies of a similar metallicity. The origin 
of this difference should be understood before using the 
MB absorbers as local counterparts of high-z DLA sys- 
tems. 

The nitrogen deficiency may be linked to the synthesis 
process. It is already known that nitrogen is produced 
by low and intermediate mass stars (LIMS; M/Mq < 8), 
while a-elements are produced by massive stars. For DLA 
systems with low nitrogen (hereafter, low nitrogen DLAs: 
LNDL As), there have been two possible scenarios pre- 
sented (jHenrv fc Prochaskall2007l and references therein): 
i.e., (i) they are in early productio n stages of n i troge n 
released from LIMS (delay scenario, iPettini et al.ll2002D . 
and (ii) they have an intrinsically flattened or truncated 
initial m ass function (IMF) w ith fewer LIMS (reduction 
sc enario. iProchaska et al.ll2002f) . However, as summarized 
in lHenrv &: Prochaskal ( 20071 ). neither scenario can explain 
the nitrogen deficiency of LNDLAs perfectly; the former 
cannot reproduce a possible bimodality of [N/g] distri- 
butio n of DLAs (jProchaska et al.l 120021 : ICenturion et al.l 
|2003|) . while the latter pro duces too much iron compare d 
with the observed amount (jLanfranchi fc Matteuccill2003f) . 

Thus, we have not yet identified the origin of the nitro- 
gen deficiency. Nonetheless, it is likely that the production 
history of nitrogen is not the same between the SMC and 
the MB. The IMF of a stellar association NGC 602 in 
the SMC (a singl e power-law with a sl ope of F ~— 1.2 for 
M/Mq = 1 - 45: ISchmalzl et al.ll2008[ ) is quite similar to 
the IMF in t he solar neighb orhood (F —1.35 for M/Mq 
= 0.4 - 10; lSalpetejlT955l) . This means that the abun- 
dance patterns of the MW and the SMC are expected to 
be similar, while they can be different from that of the MB, 
whose IMF is not necessarily the same. The MB probably 
has a local star forming history independent of the nearby 
star forming regions (i.e., the LMC and the SMC), which 
is consistent with the observed results that the young stars 
in the MB are not old enough to have escaped from the 
SMC based on their peculiar motions. 

If the MB absorbers are indeed local counterparts of 
high-z DLAs, our results would have interesting implica- 
tions: (i) DLA systems along our sight-lines to the back- 
ground quasars do not necessarily represent global prop- 
erties of absorbing structures, and (ii) the large scatter 
of metallicity and [N/a] values could be due to the inter- 
nal gradient of each DLA absorber, and they are different 
from place to place along different sight-lines go through. 
For example, we would under-estimate the global metal- 
licity (and nitrogen abundance) of the Magellanic Clouds 



if our line of sight went through only the MB. Moreover, 
there might be metallicity gradient in the MB itself; higher 
in the SMC wing than in the remaining parts of the MB 
(Lchncr et al. 2008 and r eferences th e rein). As for DLAs, 
Ellison et al.l (|2005l ) and IChen et all (|2005l ) already sug- 
gested that there could be a metallicity gradi ent as a func- 
tion o f a distance from the g alactic center. lEllison et al.l 
(|2005f) and lWolfe et~an (|2003[ ) also proposed that the DLA 
region is sampling a lower metallicity than the star-forming 
regions. Thus, high-z DLA systems also may have com- 
plex internal structures like the Magellanic Clouds and 
their neighborhood. 
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Table 1 
HST/STIS Observation Log 



QSO 


RA 


Dec 




mv 


Date 


grating 


A-covcragc 




HST Datasct ID 




(h;m;s) 


(d:m:s) 




(mag) 


(yyyy mm dd) 




(A) 


(sec) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


PKS 0312-770 


03:11:55.4 


-76:51:50.8 


0.2230 


16.1 


2001 05 12 


E140M 


1150 - 1729 


12629 


O65T01010 












2001 03 07 


E140M 


1150 - 1729 


12629 


O65T02010 












2001 10 16 


E140M 


1150 - 1729 


12650 


O65T13010 












2001 10 14 


E230M 


2132 - 2984 


6060 


O65T14010 



Table 2 

Observed and Modeled Rest-Frame Equivalent Widths 



Ti'ansition 

(1) 


Wrest (spec)" 

(mA)<^ 
(2) 


VKrest(mod)'' 
(mA)<^ 
(3) 


Lya 


< 14.73 A 


8.17 A 


Mg I A2853 


496 ± 51 


685 


C I A1329 


18 ± 4 


83 " 


O I A1302 


520 ± 7 


496 


N I A1200b 


< 311 ± 12 


278 


N I A1200a 


201 ± 20 


235 


N I A1201 


169 ± 13 


179 


Mg II A2796 


1197 ± 42 


1108 


Mg II A2803 


904 ± 43 


1063 


Mn II A1199 


< 311 ± 12 


279 


Mil II A2594 


< 49 ^ 


134 


Fe II A1261 


< 694 ± 7 


534 


Fe II A1608 


463 ± 21 


410 


Fe II A2344 


828 ± 29 


765 


Fe II A2374 


539 ± 42 


546 


Fc II A2383 


806 ± 56 


876 


Fe II A2587 


835 ± 44 


780 


Fe II A2600 


1059 ± 58 


938 


Si II A1190 


490 ± 16 


412 


Si II A1193 


495 ± 14 


444 


Si II A1260 


< 694 ± 7 


533 


Si II A1304 


451 ± 8 


423 


Si II A1527 


563 ± 9 


538 


Ni II A1317 


26 ± 4 


98 


S II A1251 


40 ± 8 


42 


S II A1254 


120 ± 7 


78 


S II A1260 


< 426 ± 7 


110 


C II A1335 


< 841 ± 8 


521 


Si III A1207 


< 999 ± 13 


300 


Si IV A1394 


78 ± 7 


101 


Si IV A1403 


73 ± 8 


57 


C IV A1548 


94 ± 12 


95 


C IV A1551 


67 ± 14 


51 


N V A1239 


< 12 ^ 


2.3 


N V A1243 


< 10 ^ 


1.1 


Ca II A3935 


142 ± 23 " 


990 



^rest-frame equivalent width or detection limit. 

^rest-frame equivalent width of best model (see 
§ 4). 

■^except for Lya for which the unit is Angstroms, 
'^blending with other lines. 

■^this would be 46.8 mA if we assume most likely 
radiation field, i.e., including extra radiation from 
the MW and the LMC with 30% escape fraction 
(see § 5). 

^not detected with 5cr. 

^detected in lower resolution sp ectrum {R ^ 
6000) taken with NTT/EMMI by fliSmoker eFaTI 



Table 3 

Best Fit Model for the MB Absorber toward PKS 0312-770 



Transition'^ 




Av" 






log Nf 


6" 


log(Z/Z0)'' 


logC/' 




riff 


sizc^ 


N deficiency"^ 


log(Si/H) 


log(N/Si) 




(km s"^) 


(km s"^) 




(K) 


(cm-^) 


(km s"^) 






(K) 


(cm-3) 


(kpc) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


H I" 


175.8 


-34.1 


0.087 


22 


19.59 


6.4 


-0.7 


-5.7 


583 


0.44 


0.029 


-1.0 


-5.2 


-0.6 


H I 


207.3 


-2.6 


0.049 


29 


19.61 


8.8 


-0.5 


-5.8 


287 


0.56 


0.025 


-0.7 


-5.0 


-0.3 


H I 


236.1 


26.2 


0.027 


46 


19.68 


11.9 


-0.7 


-5.7 


501 


0.44 


0.035 


-1.4 


-5.2 


-1.0 


O 1° 


160.8 


-49.1 






(15.00) 


(10.0) 


-1.0 


-5.1 


4450 


0.11 


0.065 


-0.6 


-5.5 


-0.2 


SiW 


216.8 


6.9 






13.16 


22.1 


-0.6 


-2.6 


13300 


0.00035 


12 


0.0 


-5.1 


+0.4 



*Name of transition that is optimized for the Cloudy model. 
'^Heliocentric velocity. 

^Relative velocity from the system center whose heliocentric velocity is 209.9 km s"""". 
'^Absorption optical depth at the line center of H I 21-cm, if column (1) is H I. 
^Spin temperature measured from H I 21-cm emission line, if column (1) is H I. 
^Best Voigt profile fit value of column density of the optimized line in column (1). 
^Best Voigt profile fit value of Doppler parameter of the optimized line in column (1). 
^Best model parameter of metallicity. 
'Best model parameter of ionization parameter. 
'Gas temperature from the best model. 

^Tolal hydrogen volume density per cubic centimeter. This depends very strongly on the radiation field's absolute strength (see § 5.5). 
'Thickness of the absorber assuming a plane-parallel structure. This depends very strongly on the radiation field's absolute strength (see § 5.5). 
"^Nitrogen Deficiency, compared to the Solar abundance pattern. 

"We adopt the same parameters as for the 236 km s~^ cloud, as an example of acceptable model. 

'^This is an example of the best fit models that have acceptable range of parameters, log{Z/ZQ) > —1.0 and \ogU > —6.0. 
^We adopt the typical metallicity of the SMC, as an example of acceptable models. 
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Table 4 

Summary of Incident Radiation Models 



Model 


EBR 


MW 


LMC 


logn.^(EBR)'^ 


logn^(MW+LMC)'' 


r 












(cm-3) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


1 


Y 






-6.06 






2 


Y 


Y 




-6.06 


-6.22 


0.69 


3 


Y 


Y 


15% 


-6.06 


-4.84 


16.6 


4 


Y 


Y 


30% 


-6.06 


-4.54 


33.1 



'^Volum e density of hydrogen-ionizing photons from the extragalactic background 
radiation (|Haardt fc Madaulll996l . l200lh . 

'^Volume densi ty of hydrogen-ionizing photons from the MW and the LMC 
l^eEaD [20053) ■ 

"^Ratio of the ionizing photon number densities from the MW and/or the LMC, 
to that from the EBR, n^{MW+LMC)/n^(EBR). 



Table 5 

Model Constraints for the Magellanic Bridge Systems 



Cloud 


parameter 


constraint 


line 




condition 


(1) 


(2) 


(3) 


(4) 




(5) 


H I {Vq=207 km s^^) 


log (7 


> 


-6.0 


C I 


to 


avoid over-production 






< 


-4.8 


C I 


to 


avoid under-production 




log 


.z 


> 


-0.5 


S 11 


to 


avoid under-production 








< 


-0.5° 


Ni 11 


to 


avoid over-production 


H I (Vq=236 km s"1) 


log 


lU 


> 


-6.1 


C I 


to 


avoid over-production 








< 


-5.0'' 


Si II 


to 


avoid under-production 




log 


\Z 


> 


-0.8 


Si II 


to 


avoid under-production 








< 


-0.6° 


Ni II 


to 


avoid over-production 


H I (^0=176 km s-^) 


log 


lU 


> 


-6.0" 


C I 


to 


avoid over-production 






< 


-5.0" 


Si II, O I 


to 


avoid under-production 




log 


\Z 


> 


-1.0 


Si II 


to 


avoid under-production 








< 


-0.7 


Mg I, Ni II, C I 


to 


avoid over-production 


Si IV (V0=217 km s"^) 


log 


■.V 


> 


-2.7 


C IV 


to 


avoid under-production 






< 


-2.4 


C IV 


to 


avoid over-production 




log 


.z 


> 


-4.0"* 


H I 21-cm 


to 


avoid over-production 


O I (^0=161 km s-^) 


log 




> 


-6.0 


Mgl 


to 


avoid over-production 






< 


-3.5 


Si IV 


to 


avoid over-production 




log 


\Z 


> 


-1.0"* 


H I 21-cm 


to 


avoid over-production 



^Because nickel is strongly depleted on to the dust, this upper limit on the metallicity can be softened 
once the depletion is considered. 

^This constraint is valid if we assume \og{Z/ Zq) — —0.7. 

*^This constraint is valid if we assume —1.0 < \og{Z/ZQ) < —0.7. 

■^No upper limit can be placed because the corresponding H I 21-cm line is not detected with > Icr 
detection limit. 
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Table 6 

Comparison of Abundance Patterns" 



\og{elernent / H) 





C 


N 


O 


Si 


S 


ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


The Sun 


-3.61±0.04 


-4.17±0.11 


-3.31±0.05 


-4.46±0.05 


-4.80±0.05 


1 


SMC H II region 


-4.47±0.06 


-5.41±0.08 


-3.95±0.08 


-5.30±0.2 


-5.58±0.11 


2 




-4.84±0.04 


-5.54±0.12 


-3.98±0.08 




-5.51±0.14 


3 


SMC star (AV 304) 


-5.15 


-5.11±0.17 


-3.90±0.16 


-5.27±0.02 




4 


Bridge star (DI 1162) 


-4.94±0.33 


-5.30 


-4.00±0.13 


-5.19±0.06 




4 


Bridge star (DGIK 975) 


-5.32±0.25 


-5.26 


-3.96±0.21 


-5.78±0.15 




4 


Bridge gas {Vq — 176 km s^^) 


-4.61 to -4.31 


-6.37 to -5.37 


-4.31 to -4.01 


-5.46 to -5.16 


-5.74 to -5.44 


5 


Bridge gas (Vq = 207 km s^^) 


-4.11 


-5.37 to -5.17 


-3.81 


-4.96 


-5.24 


5 


Bridge gas (Vq — 236 km s^^) 


-4.41 to -4.21 


-6.47 to -5.87 


-4.11 to -3.91 


-5.26 to -5.06 


-5.54 to -5.34 


5 


Bridge gas (toward DI 1388) 




-5.29±0.11 


-4.30+0.13-0.11 


-5. 45+0. 14-0. 12'^ 




6 



^AU abundance patterns are measured from the observed spectra directly, except for our results for which we estimate them based on the 
photoionization model. 

^1: ILodderd ( |i>00j|) . 2: IKurt et all UOO^) , 3: IDufouH l IlOM . 4: IHolleston et al.| (119^)911 . 5: this paper, 6: ILehner et aTl diiOOBIl 

■^converted from log(0/H), assuming solar abundance pattern. 
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Galactic Longitude 



Fig. 1.— Location of our target, PKS 0312-770, as well as stars in the MB (DGIK 97 5 and DI 1162, and DI 1388) and the SMC (AV 304), 
superimposed on the H I peak brightness temperature (Tg) from iPutman et al.l II1998I) . Contours mark the 0.1 K (5 it), 0.8, 2, 8, 16, 32, 
64, and 128 K peak brightness. Light and deep gray contours correspond to Tb >2 K and >8 K, respectively. Light gray contour roughly 
corresponds to the minimum threshold DLA column density, Nuj = lO'^" '^ cm~^ if the total line width of ionHl 21-cm absorption is similar 
to that toward PKS 0312-770, dv ~ 50 km s'^ (see Figure 5). 
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1302.6 1302.8 1303 1303.2 1303.4 



Observed Wavelength (A) 



Fig. 3. — M odeled spectrum with three O I A1302 components whose Doppler parameters are b = 5.8, 7.6, and 6.3 km s~^, from left to right 
l|Kobulmckv fc Dickev 1999). Positions of the lines are marked with ticks. Even if we synthesize spectra varying the column density of each 
component in the range log Af (O l/[cm~^]) = 15 — 19 (i.e., a total column density of all three components is in the range log Af (O l/[cm~^]) 
~15.5 — 19.5), both sides of the observed spectrum shape cannot be reproduced. We test a ridiculously wide range of column density by using 
same values for three O I components without any justifications, because the purpose of this test is only examine whether the given Doppler 
parameters could reproduce both sides of the observed absorption profile. 




1302.6 1302.8 1303 1303.2 1303.4 



Observed Wavelength (A) 



Fig. 4. — Same as Figure 3, but the three O I A1302 components have Doppler parameters of 6 = 6.4, 7.6, and 11.9 km s ^ from our own 
fitting trials. A model with logA^ (O l/[cm-2]) of 16.0 gives an acceptable fit at the high-velocity side of the observed spectrum. 
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Fig. 5. — Modeled spectr um of the H I 21-cm absorpt ion lines, overlaid with the radio spectrum of PK S 0312—770 taken with the Australia 
Telescope Compact Array l|Kobulnickv &: Dickeviri999t ). In addition to the original three components ijKobulnickv Dickeviri999l ). we add 
another component at At; = —49 km s~^ (i-e., Yq = 160.8 km s~^), required by the modeling of the metal lines. This additional component 
(whose central optical depth is t = 0.022) is not detected at more than a \a level in the observed spectrum whose S/N ratio is about 40 per 
pixel (which corresponds to the detection limit of t = 0.025 at the line center). 
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Fig. 6. — Strength of the ionizing radiation field as a function of distance along our sight-line toward PKS 0312—770. The solid curve 
includes the radiation from the LMC as well as the radiation from the MW, while the dotted curve includes only the latter. The black 
curves are for the direction to PKS 0312—770, while the red curves are for the direction straight toward the LMC. For example, at a radius 
of 50 kpc, the ionization strength decreases from (i) LMC + MW toward LMC, (ii) LMC + MW toward PKS 0312-770, (iii) MW toward 
PKS 0312-770, to (iv) MW toward LMC. The LMC model assumes that 30 % of the ionizing radiation escapes from the LMC, which we 
consider to be an upper limit. 
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Fig. 7. — Spectrum of four different radiation fields that we applie d in our photoioni z ation calcul ations with Cloudy. The smooth, dashed, 
red line denotes the extragalactic background radiation (EBR) from lHaardt &: Madaul II19961 . 120011) . The three solid blue lines represent the 
radiation fields from the MW, MW + LMC with 15% escape fraction, and MW + 30% escape fraction, from bottom to top. These are 
normalized assuming a distance of D = 50 kpc from the center of the MW to the PKS 0312—770 sightline cloud. Ionization edges of several 
important transitions arc indicated at the top of the plot. 
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Fig. 8. — Detected transitions, and tiiose that provide limiting constraints, shown in velocity space for the absorption system in the 
Magellanic Bridge toward the quasar, PKS 0312-770. The velocity range is ±300 km s~^ for all transitions except for Lya whose range is 
±3000 km s~^. Velocity = km s~^ corresponds to the system center of the Magellanic Bridge absorber whose heliocentric velocity is Vq 
= 209.9 km s''^ . The data are from a HST/STIS spectrum. The error spectrum is also indicated as a solid histogram just above the dotted 
line crossing each plot at zero flux. An example of the best model fit to the observed spectrum (using model parameters listed in Table 3) is 
superimposed on the data as a solid (red) curve. Three unshaded regions denote the locations of the Milky Way, the MB, and a high-velocity 
cloud (from left to right), respectively. The positions of absorption components are marked with ticks. Above the spectrum in each panel, 
three ticks in the bottom line are H I optimized components, while those in the top line are O I or Si IV optimized components. 
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Fig. 8. — Continued. 
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Fig. 9. — Comparison of the model spectra with nitrogen deficiency (solid red line) and without nitrogen deficiency (dashed blue line). 
Obviously, the model with solar abundance pattern does not reproduce the observed N I absorption lines. There is no large difli'erences in 
the N V doublet at this ionization parameter because of their very weak strengths. Unshaded parts are the regions of the absorber in the 
Magellanic Bridge. 
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Fig. 10. — Comparison of the model spectra using different incident radiation fields: (a) only extra-galactic background radiation (EGR) 
(red), (b) EGR plus radiation from the MW at a distance of D = 50 kpc from the MW (green), (c) EGB, MW radiation, plus radiation from 
the LMC with a 15 % escape ratio as a minimum flux model (thick blue), and (d) EGB, the MW radiation, plus radiation from the LMC 
with a 30 % escape ratio as a maximum flux model (thin blue). Unshaded parts are the regions of the absorber in the Magellanic Bridge. 
The color version of this figure is available in the electronic edition of the Journal. 
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Fig. 11. — log(N/Si) versus log(Si/H) of three absorption components in the Magellanic Bridge toward P KS 03 12—770, compared to 
those of damped Lyman alpha systems (dots; Henry & Prochaska 2007), blue compact galaxies (b lue cross: INava et al. 2006) in which 
oxygen abundance was converted to silicon abundance using the solar abundance ratio ( H olwegerii200ll : lAUende Prieto et al.i.2001i) , and other 
SMC/MB objects from Table 6 (error bars). Error bars shown in the bottom left are average uncertainties of the DLA systems in each 
direction. The dotted error bar denotes the la error in our determination of the abundance ratio for the H I component at Vq = 236 km s~^, 
if we reduce the nitrogen abundance by —1.0 dex instead of —1.4 dex that is still consistent with the observed spectrum. The four possible 
evolutionary vectors are also shown in the bottom right. We do not plot components at V@ = 161 km s~^ and 217 km s~^, because we cannot 
place any meaningful constraints on their abundance patterns. 
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Fig. 12. — Close up velocity plot around three detected N I absorption lines within ±100 km s~^ of the system center. Solid lines represent 
a model in which we decrease nitrogen abundance of the H I component at Vq = 237 km s~^ by 1.4 dex compare to the other elements (i.e., 
our best model in Table 3), while dashed lines are for the same model but with the nitrogen abundance reduced by 1.0 dex. If the continuum 
fitting around N I 1200a is incorrect, this latter model would reproduce absorption features of the other two N I lines much better. Moreover, 
this weaker nitrogen deficiency is more reasonable if we compare to those of the other two H I components at Vq = 176 and 207 km s~^. 



